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EXPERIMENTAL INVESTIGATION OF HEAT -TRANSFER 
AND FLUID-FRICTION CHARACTERISTICS OF 
WHITE FUMING NITRIC ACID 
By Biruce A. Reese and Robert W. Graham 

SUMMARY 




As part of the general rocket-research program of the National 
Advisory Committee for Aeronautics-, e3iperiments have been conducted to 
determine the heat-transfer and fluid -friction characteristics of white 
fuming nitric acid over the following range of conditions : 


Heat-flux density, Btu/(sq in.)(sec) 0.13 - 1.^ 

Pressure, Ib/sq in. abs 6k - I 65 

Acid inlet temperature, °F 50 - 137 

Reynolds number 55,000 - 220,000 


In the forced-convection regime (without surface boiling) the 
Fanning friction coefficient with heat transfer f can be satisfac- 

torily related to the Isothermal friction coefficient f^go measured 
at the same Reynolds number by the following eq-uation; 


^Iso _ / ^^bulk \*^ ’ ^ 

^q X^Vall/ 

where q is the viscosity of the white fuming nitric acid. 

The heat-transfer results are satisfactorily correlated by the 
following equation; 


J = 0.023 

where j is the -Colburn factor, N|^g is Reynolds number, and all 
physical properties are evaluated at the bulk temperature. 

Some preliminary results are presented for the regime of forced 
convection with nucleate boiling. 
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INTRODUCTION 


In recent years a number of investigations have been conducted on 
heat transfer to liquids where large temperature gradients occurred at 
the heated surface. (Some of these investigations are reported in 
refs. 1 to 16.) In the design and operation of- regenerative ly cooled . 
rocket motors it is desirable to have an accurate knowledge of the heat- 
transfer and pres sure -drop characteristics at the liquid coolant. The 
correlations based on convective heat-transfer measurements at low heat- 
flux densities have not been considered adequate for the prediction of 
the convective heat transfer at the high heat -flux densities encountered 
in regeneratively cooled rocket motors. The heat-transfer character- 
istics of a few of the rocket propellants employed as coolants have been 
investigated and reported in the literature (refs. 1 to 3), 

As part of- its general rocket-research program the National Advisory 
Committee for Aeronautics initiated an experimental- investigation at 
Purdue University Rocket Laboratory on the behavior of white fuming 
nitric acid (commercial anhydrous nitric acid, hereafter designated as 
WFNA) as a coolant and the results obtained are reported herein. This 
investigation was carried out under the sponsorship and with the financial 
assistance of the NACA. An investigation of the heat-transfer and 
pressure-drop characteristics of WFNA was particularly desirable because’ 
the WFNA decomposes even at room temperatures. Consequently it was not 
at all certain that the heat-transfer correlation equations derived for 
pure liquids at low heat-flux densities could be employed for calculating 
the heat-transfer coefficients for WFNA at high heat-flux densities 
because of the possibility of the composition of the WFNA changing with 
heat addition; in decomposing the WFNA yields oxides of nitrogen and 
water. From the more fundamental point of view of studying heat transfer 
per se, WFNA is an unsatisfactory liquid which introduces appreciable 
experimental difficulties because of its corrosivity and uncertainties 
regarding its physical properties, 

A 0.539-inch- inside-diameter tube 24 inches long fabricated from 
Haynes-Stellite Alloy 25 was utilized as the test section, which was 
employed as a resistance element heated by alternating current and cooled 
by the WFNA flowing throu^ it. Experimental data were taken over a wide 
range of operating conditions: Resmolds numbers from 55,000 to 220,000, 

heat-flux densities ranging from O.13 to 1.4 Btu/(sq in.)(sec), system 
pressures from 64 to 165 Ib/sq in. abs, and inlet WFNA temperatures from 
51° to 137° R- Most of the heat-transfer data were in the forced- 
convection (without surface boiling) regime, but some preliminary inves- 
tigations were made In the forced-convection region where the surface 
temperature of the test section exceeded the ..indicated boiling point 
(surface bailing) of the WFNA. Accurate data on heat transfer with sirr- 
face boiling are of considerable -interest since investigations with water 
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In this regime indicate that heat-flux densities up to 8 Btu/(sq in.)(sec 
can he absorbed by the coolant. Heat-transfer rates of this order of 
magnitude are now being encountered in the research being conducted at 
the Purdue University Rocket Laboratory concerned with the operation of 
rocket motors at abnonnally . high combustion pressures (ref. 17 ). 

Acknowledgment is due Dr. M. J. Zucrow, project director, for his 
continuous assistance throu^out the project and for his help in the 
preparation of this report. The authors would also like to express 
their gratitude for the wholehearted cooperation they have received from 
the NACA staff, both in Washington, D. C., and at the Lewis Flight 
Propulsion Laboratory, Cleveland, Ohio. 


SYMBOLS 


A eirea, sq ft, sq in. 

C constant 

Cp specific heat at constant pressure, Btu/(lb)(°F) 

D diameter, ft, in. 

f Fanning friction factor 

G wei^t rate of flow, Ib/hr, Ib/sec 

g gravity constant, 32.17^ ft/sec^ 

h surface coefficient of heat transfer, Btu/(sq ft)(hr)(°F), 

Btu/(sq in. )( sec)(°F) 

J Colburn factor 

k thermal conductivity, Btu/(ft)(hr)(°F), Btu/( in. )(sec)(°F) 

L length, ft 

M wei^t flow rate per unit area, lb/(hr)(sq ft) 

Nusselt number, ho/k 

Wpj. Prandtl number, CpM-/^ 

Npg Reynolds number, pDV/p. 
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Ap frictional pressure drop or pressure drop across orifice, 

Ib/sq. in. 

P pressure, Ib/sq in. abs 

q heat transferred per unit time, Btu/hr, Btu/sec 

Pi radius of inside of tube, ft 

'Bq radius of outside of tube, ft 

t bulk or mixing-cup ten5>erature, °F 

to outside siirface temperature of test section, °F 

tg surface or inside tube-wall temperature, °F 

'^'sat saturation temperature, °F 

Atx excess tengjerature, tg - tg^^, °F 

V velocity, ft/sec 

X distance from inlet 

y specific weight, Ib/ft^ 

p viscosity, lb/(ft)(sec) 

p density 

Subscripts: 

bulk properties evaluated at bulk temperature 

iso isothermal 

m metal 

mean mean value 

q nonisothermal 

s inside surface or wall of tube 

wall measured at inside surface of wall 
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1 

inlet 

2 

outlet 

Superscripts: 

m 

constant 

n 

constant 


APPAEATUS AND INgEEUMENTATION 


The design and arrangement of the apparatus and instrumentation 
employed in the investigation are discussed in detail in reference 18. 
Only a brief description is, therefore, presented here. 

Figure 1 is a schematic diagram of the electrically heated test 
section and the associated equipment. Referring to figure 1 the flow 
throu^ the apparatus is as follows. The pump A circulates the acid 
in the closed system. The flow was regulated by the remotely controlled 
air-operated Annin throttle valve C and the hand-operated bypass valve B. 
The flow was metered by the calibrated sharp-edge orifice D with the 
differential pressure being determined by means of a Wiancko differential 
pressure transducer and a hand-balanced potentiometer. . The small tank E 
was Installed for reducing pressure and flow fluctuations and for con- 
venience in pressurizing the flow circuit. 

Two Annin throttle valves F and J were Installed at the upstream 
and downstream ends of the test section as a safety feature; controls 
were Incorporated for closing the valves rapidly in the event of a test- 
section rupture. The two mixing chambers located upstream and downstream 
from the test sections con5>rlse three concentric tubes arranged so that 
the acid makes three passes longitudinally through the mixipg chamber 
before its temperature is measured. This design of the mixing chamber 
assures that the thermocouple reads the bulk (mixing cup) temperature. 

To obtain a fully developed isothermal boundary layer and a iiniform 
tiirbulent velocity profile for the WFNA entering the test section a 
starting length of tube equal to approximately 25 tube diameters was 
inserted between the upstream mixing chamber and the test section. A 
strai^t length of tube equal to approximately 10 tube diameters is 
employed between the test section and the downstream mixing chamber to 
insiire smooth flow leaving the test section. 

The test section H is a Haynes-Stellite Alloy 25 tube having an 
outside diameter of 5/8 inch, a wall thickness of O.OU3 inch, and a total 
length of 24 inches. Figure 2 presents the details of the Installation 
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of the test section In the WMA flow circuit. The outside wall ten5)era- 
tures of the test section were measured at l 4 locations (see fig. 2(c)) 
with Chrome 1 -Alumel themocouples (30 gage wire) and an automatic self- 
balancing indicating- type Brown potentiometer. The pressure drop over 
the length of the test section was measured by utilizing a Wlancko dif- 
ferential pressure transducer. The test section was insulated thermally 
by a l/ 2 -lnch layer of Sauereisen cement and 12 inches of glass wool with 
two concentric aluminum- foil radiation shields spaced at radii of 2 and 
6 inches approximately. The glass wool and aluminum also surrounded the 
tubes which connected the mixing chambers to the test section. 

The test section was heated with alternating current, the tube 
acting as a resistance element, and the WMA cooled the test section by 
forced convection. Power was supplied to the test section from a 
2ltO-volt alternating-current single-phase generator and a multitap trans- 
former having a maximum power rating of 100 kilowatts . Because of the 
maximum current limitation of the raultltap transformer (2,500 amps) anri 
the resistance of the test section the maximum utilizable power for 
heating was 67 kilowatts. 

The heat added to the WPNA after flowing through the test section 
was removed in the double-tube aluminum heat exchanger K, which enabled 
one to maintain a constant tenperature for the WFNA entering the test 
section. Valves M, N, and P were provided for draining the system and 
flushing it with either air or water. The system is filled with WPM 
by forcing the WFNA out of the storage tank Q with low-pressure nitrogen 
gas. 


The WFNA flow system and instinimentatlon described in reference 18 
performed satisfactorily in preliminary trials made with water as the 
test fluid. After some exploratory experiments with WFNA it was neces- 
sary to make the following changes: 

(1) The thermocouples in the mixing chambers were changed from 
Chromel-Alumel to copper-constant an thermocouples (30 gage)j the copper- 
constantan couples were found to be more suitable for the range of WFNA 
temperatures encountered. 

(2) The thermocouples were silver-soldered to the outside siu*face 
of the test section instead of using the condenser discharge method for 
attaching them; the condenser discharge method caused the thermocouples 
to become so brittle that they broke quite readily. The silver-solder 
method consisted of drilling a hole O.OO5 inch deep with a number 40 drill, 
filling the depression with silver solder, melting the silver solder, and 
then quickly inserting the thermocouple in the molten solder. This method 
for attaching the thermocouples produced a strang junction at the surface 
of the tube and had no discernible effect upon the properties of the 
thermocouple wire. 
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(3) The annular-ring method for measuring the static pressure at 
the entrance and exit of the test section had to be abandoned because of 
electric current being conducted by the WFWA across the small gap between 
the flanges which caused the weld that attached the Van Stone ring to 
the tube to burn through, thus producing leakage of WFIfA. The pressure 
taps employed in all of the experiments consisted of l/l 6 -inch holes 
drilled throu^ the flanges of the mixing chambers that were attached to 
the test section. Figure 2 (a) shows the position of the drilled holes 
and Indicates by means of dotted lines the discarded method for measuring 
the static pressure. The drilled holes are located approximately 1/2 inch 
from the ends of the test section, but the corrections for the unheated 
portion of the connecting tubes are relatively simple to make. 

( 4 ) After the method for measuring static pressure was changed 
(see (3) above) the life of the test section was increased. It was 
found, however, that electric current was still being conducted through 
the WFNA from the test section to the mixing chamber and burning through 
the weld that attached the Van Stone ring. To reduce the flow of electric 
current the mixing chambers were Insulated electrically from the rest of 
the WFNA flow circuit. Electrical insulation of the mixing chambers from 
the rest of the circuit reduced the rate of burning at the ends of the 
test section by arcing, but did not eliminate it entirely. 

(5) The Statham differential pressure transducers, which functioned 
satisfactorily for a long period of time over a wide range of conditions 
when water was used as the coolant, failed to operate aJter only short 
exposure to the WFNA probably because of corrosion. Replacement Statham 
units also failed after short exposure to the WFNA. Since Wiancko dif- 
ferential pressure transducers had been used successfully at the Purdue 
University Rocket Laboratory for WFNA flow and pressure measurements, the 
Statham transducers were replaced by Wiancko transducers . The Wiancko 
transducers performed satisfactorily in the WFNA circuit, the only dif- 
ficulty encotintered being slight change in the zero (no differential pres- 
sure) reading with time, presstire, and temperature; consequently it was 
necessary to read the zeros at frequent intervals. 


PRELIMINARY TESTS AND PROCEDURE 


A nvimber of preliminary investigations were conducted to check the 
Instrumentation . 


Measurement of Weight Flow 

The Statham differential pressure transducers were calibrated against 
a mercury manometer. Employing the Statham unit as the differential 
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pressure gage, the sharp-edge orifice in the WFM flow circuit was 
calibrated' in position. Water from the cooling-system supply tank was 
pumped through the WFNA flow circuit and weighed as it discharged from 
that circuit. 


Measurement of Bulk Temperature 

The copper-constantan thermocouples employed for the measurement 
of the bulk temperature of the WFNA were calibrated in a furnace against 
a standardized platinum - platinum- rhodium thermocouple. The deviations 
from the standard were negligible. _ 


Measurement of Surface Temperature of Test Section 

The temperatures of the outside surface of the test section were 
read on a 20-point-indicatlng Brown automatic self-balancing potentiom- 
eter. When the tube was heated with alternating-current power it was 
found that the indicating dial of the Brown potentiometer slowly 
oscillated (approximately l/2 cycle per second) when connected to those 
thermocouples attached at the ends of the test section; the amplitude 
of the oscillations increased with increasing temperature to approxi- 
mately ±20° F at 500° F. The thermocouples near the center of the test 
section, however, always gave steady readings. The temperature readings 
of the thermocouples located at the end of the test section when measured 
with a hand-balanced potentiometer were always steady. Oscilloscope 
measurements of these same thermocouples showed both a 60- and 900-cycle 
alternating -current pickup from the thermocouples. It was not possible 
to eliminate the oscillations of the Brown potentiometer dial but it was 
fovind experimentally that the average of the maximum and minimum readings 
of the dial gave the correct values of aurface temperature. 

In an effort to locate the source of the oscillations the magnetic 
field around the test section was investigated by placing a loop of wire 
connected to a sensitive galvanometer at a number of positions. No 
field was detected that had sufficient strength to Influence the tempera- 
ture readings . 


Heat Losses 

The heat loss from the insulated test section was obtained for two 
different average wall tenperatures by supplying power to the test 
section without liquid in the system. After equilibri;im conditions had 
been maintained for 1 hour the power input and all test -section outer 
surface temperatures were recorded. At an average surface temperature 
of 2kO° F the power loss was 5 watts and at an average surface temperature 
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of 450 ° F the power loss was 11.3 watts. The corresponding heat loss 
was always less than 0.1 percent of the power supplied to the test sec- 
tion during any heat-transfer experiment and was, therefore, neglected 
in calculating heat-transfer coefficients. 


Heat-Transfer Test With Water 

Heat-transfer tests with water were conducted to test the apparatus 
and the experimental procedure. The results, which covered a Reynolds 
number range of 4 - 9,500 to l 66,000 with heat-flux densities up to 
0.76 Btu/(sq in.)(sec), were satisfactorily correlated by the equation 
of Kaiifitian and Isely (ref* 4) : 

%u = 0.0168 (NRe)°‘^^ (Npr)°'^ 


Heat-Transfer Tests With WFNA 

Before conducting the heat-transfer experiments with WFNA, the modi- 
fications to the apparatus discussed in the section entitled "Apparatus 
and Instrumentation" were made. The detailed steps of the experimental 
procedure were as follows: (l) The cooling water was circulated through 

the aluminum heat exchanger; (2) the acid system was filled and pres- 
siirized to the desired value; (3) the zero readings were taken for the 
Wiancko units; (4) the acid pump was started and the acid flow was 
regulated to a predetermined value; and (5) the motor generator set was 
started and electric power was applied to tlie test section. 

After thermal equillbri\im was established, readings were taken of 
the outside s\irface ten 5 >eratures of the test section, the bulk tempera- 
tures of the acid entering and leaving the test section, the electrical 
power supplied to the test section, the inlet pressure to the test sec- 
tion, the acid flow rate, and the pressxire drop across the length of the 
test section. 

In most of the experiments the electrical power supplied to the test 
section was varied while the WFNA flow rate and inlet temperature were 
maintained relatively constant. At the completion of each series of 
such trials the electrical power and the acid pump were shut down and 
the zeros of the Wiancko units were read. For the tests with the higher 
acid inlet temperatures the flow of the cooling water in the aluminxim 
heat . exchanger was reduced until the acid temperature weis raised to the 
desired value by the electrical power supplied to the test section. The 
cooling-water flow rate was then adjusted until thermal equilibrium was 
established at the desired WFNA inlet temperature. 
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METHOD OF CALCULATION 


The calculations In this investigation vere carried out on IBM 
machines at the Purdue Statistical Laboratory. The method of calcu- 
lation, the important equations, the physical -property data, and an 
estimate of the accuracy of the results. are given below. 


Temperatures 

Since the following heat-transfer coefficients were desired, 

(a) average over-all value for the test section and (b) the local values 
at seven positions along the test section, it was necessary to compute 
( 1) the average temperature of the outside surface of the test section 
and (2) the average temperature for the two .thermocouples employed for 
measuring the local value of the surface tenqjerature at each of the seven 
positions. The average temperature of the outside surface of the test 
section, denoted by tQ, was dete2miined by integrating numerically the 

measured ten?)erature along its surface as a finaction of the test-section 
length. The temperatures used in the integration were the average of 
the two thermocouples located l8o° apart at each of the seven stations 
along the test section. 


The average temperat\ire of . the inside surface of the test section 
denoted by tg was calculated from the following equation (ref. 19): 


^s = 




( 1 ) 


It was assumed in deriving equation ( l) that the heat was generated uni- 
formly across the tube wall and that all of the heat flow was radially 
inward. Substituting the values of = O.625 inch, % = 0.539 inch, 
and L = 2 h inches into eqiiation ( l) gives ^ 


^s ” ^o 


0.00617 ^ 


( 2 ) 


where q is the electrical power supplied to the test section in British 
thermal units per hour. 

Equations (l) and (2) may be corrected to Include the effect of the 
variation of electrical resistivity with temperature (see refs. 5 and I8), 
but since these data are not available it was not felt that the experi- 
mental work necessary to obtain these data was justified because of the 
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absence of precise data for the thermal conductivity of the metal and 
for the physical properties of WFKA. 

The heat loss from the test section was so small ( see the section 
entitled "Preliminary Tests") that the measured value of q, with no 
correction for heat loss, was utilized as the heat transferred to the 
WFNA. The values of thermal conductivity of Haynes -Stellite Alloy 25 
used in equation (2) are presented in figure 3 * The thermal conduc- 
tivity km of the alloy was determined experimentally at lk5° F. The 
slope of the curve from the measured value was extrapolated from pub- 
lished data on the thermal conductivity of several other similar Haynes - 
Stellite alloys (ref. 20). The temperature employed for determining the 
value of km from figure 3 was the average tube-wall ten^ierature 
(to + tg)/2. 


Equations were derived for the curves presenting the thermal conduc- 
tivity of Haynes-Stellite Alloy 25 as a function of temperature and all 
the physical properties of WFNA as a function of temperature. These 
equations were converted into coded instructions for use in the IBM 
machine calculations. 

The average bulk temperature t of the WFNA utilized in calculating 
the average over-all heat-transfer coefficient and the Fanning friction 
factor was assumed to be the arithmetic mean of the bulk temperatures of 
the acid entering the test section t]_ and that leaving the test 
section tp. 

In calculating the local heat-transfer coefficients the arithmetic 
average of the readings of two thermocouples l 80 ° apart was used for toJ 
the inside surface temperature tg was calculated by means of equa- 
tion ( 2 ). In calculating the bulk temperature of the acid for each of 
the seven stations it was ass\amed that the WFNA temperature Increased 
linearly from the inlet to the outlet of the test section. 


Heat-Transfer Coefficient 

The heat-transfer coefficient h was calculated from the equation 


h = 


<1 

A(t8 - t) 


( 3 ) 


where q is the electric power supplied to the tube in British thermal 
units per hour and A is the inside surface area of the test section 
in square feet. 
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Fanning Friction Coefficient 

The friction data were evaluated from the observed pressure drop 
across the length of the test section with and without heat addition. 

The momentiim pressure loss 'across the test section was Insignificant; 
thus, no attempt was made to correct the observed pressure-drop readings 
for the momentum pressure loss. Average Fanning friction coefficients 
for the test section were calculated from the observed data by employing 
the following equation: 

f = ApgyA^P 
2G^L 


Physical Properties of WFNA 

All of the physical properties of WFNA employed in calculating the 
dimensionless groups Nj^g, ^Pr> ^Nu» with the exception of the 

specific heat, were determined at Purdue University under the sponsor- 
ship of the NAGA (see ref. 2l). The physical -property data are j^resented 
in figures 1 to 4 of reference 21 as functions of temperature. The 
Prandtl modulus is presented in figure 6 of reference 21. The error in 
the values of the specific heat of WFNA is unknown. The maximum errors 
in the viscosity meas-urements are j-udged to be of the order of 5 percent 
at 300 ° F and less than 1 percent at room temperature. The error in the 
density determinations is Judged to be approximately 0.2 percent in the 
region below 100° F and approximately 1; percent at 300° F. The thermal- 
conductivity measurements are Judged to be within 5 percent (up to 122° F) 
of the correct value. The extrapolation of the latter data to 300° F is 
probably of the same order of accuracy. 


Estimated Accuracy of Calculated Results 

The factor contributing the greatest probability of error in the 
values of the heat-transfer coefficient is the determination of tlie 
inside surface temperature of the test section. The thermocouples were 
silver- soldered in a small hole approximately O.OO 5 inch deep. It is 
probable that the temperature of the silver solder which is measured by 
the thermocouple is somewhat lower than the outside temperature of the 
test section. The maximum probable error in measuring the outside sur- 
face tenperature of the test section is approximately -4 percent of the 
temperature drop through its wall. Experiments were conducted that 
justified the assumptions made in deriving equation (2), which was employed 
for calculating the temperature drop through the wall of the test section. 
There are, however, errors in the values of q and k that were substi- 
tuted into equation (2). The error in measuring the electrical power q 
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was of the order of +1 percent, hut the leakage of the electric current 
to the surroundings, which was not taken into account, is estimated to 
be as much as 2 percent. The value of kjjj for the test-section material 
is estimated to be uncertain to ±7 percent. Consequently the error in 
the calculated value for the inside surface temperature of the test sec- 
tion might be as large as l4 percent of the temperature drop through the 
tube wall. This is deemed to be the maximum probable error in the 
value tg - tg. 

The \mcertainty in the value of the heat -transfer coefficients cal- 
culated by n^ans of equation ( 3 ) is primarily a function of the error in 
the inside surface temperature of the test section, the inside surface 
area and bulk temperature being known to within +1 percent. For example, 
consider run 21.4 (table l) . An uncertainty of l4 percent of the tempera- 
ture drop through the tube wall causes an error of 7° in the calculated 
temperature difference (ts - t) . Since the measured tenrperature differ- 
ence was 120°, the maximum probable error in the heat -transfer coef- 
ficient h is ±6 percent. 

The probable error in measuring the weight flow rate of WFNA is 
±4 percent and is due to the following uncertainties: (a) Uncertainty 

caused by the scatter of 2 percent in the coefficient of discharge of 
the orifice plate, and (b) uncertainty of 2 pertent in the values of the 
density of WFNA due to changes in its composition. 

The probable error in the Fanning coefficient of friction is a 
function of error in the weight-flow measurement, the measured pressure 
drop, and the acc^lracy of the specific-weight measurement. The error in 
the pressure-drop measurement is estimated to be less than 3 percent and 
that of the specific weight is 2 percent. These uncertainties give an 
estimated maximum error of ±7 percent in the value of the Fanning friction 
coefficient . 

The uncertainties in the values of the j -modulus and are 

primarily due to the uncertainties in the physical-property data for WFNA. 
The j -modulus is a function of the following variables: 

j = f(h,G,Cp^/3^j^2/3^j^2/3^ 

Based upon the probable errors presented above for the properties 
of WFNA the uncertainty in the values of the J-modulus could be as much 
as 16 percent assuming that the composition of the WFNA was the same as 
that used in determining the physical-property data. Even larger uncer- 
tainties are possible if the congjosition of the WFNA used in the heat- 
transfer experiments differed from that used in' determining the physical 
properties. 


1 
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The Reynolds number is a function of the following variables : 


The maximum probable error in the values of N^g based upon a WMA 
of the same composition as that used in determining its viscosity is 
±5 percent. For WFNA compositions that are . different the uncertainties 
would be larger. 


RESULTS AKD DISCUSSION 
Results 


Heat -transfer data for WFNA flawing thirough a 24- inch-long Haynes - 
Stellite Alloy 25 tube ( 0.539-inch inside diameter) were obtained over 
the following range of variables; 


Heat-flux density, Btu/(sq in.) (sec) 0.13-1.4 

Pressure, Ib/sq in. abs 64 - 165 

WFNA inlet temperatixre, °F 50 - 137 

Reynolds number 55^000 - 220,000 


The results are presented in table 1. 


Heat Balance 

The heat transferred to the WFNA was determined by two methods: 

(a) By measuring its weight flow rate and its buHc-ten^jerature change 
and then applying the equation q = Gcp(t 2 - "tl); and (b) by measuring 

the electric input with a wattmeter. If heat losses are neglected these 
two measurements should be equal. Figure 4 is a plot of the heat input 
to the test section measured electrically as a function of the heat trans- 
fer based on method (a). It can be seen from figure 4 that most of the 
measurements agree within ±5 percent. The large deviations which occurred 
at small values of bulk -temperature difference are believed to be due to 
inaccuracies in measuring the bulk teniperature. The inaccuracies in 
measuring the weight rate of flow of WFNA, combined with ;incertainties in 
its specific-heat data, probably account for some of the scatter in the 
data. There was also some conduction of electrical current throu^ the 
WFNA that could not be meastired. The heat balance Indicates that the 
current conducted throu^ the WFNA was a small fraction of the total 
current supplied to the test section. 
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Correlation of Average Fanning Friction Coefficient 


Within the convective heat-transfer regime, the Fanning noniso- 
thermal^ friction coefficient fg is consistently less than the Fanning 
isothennal friction coefficient f^so measured at the same flow rate. 

For the range of average bulk temperatures listed above, the noniso- 
thennal friction coefficient is related to the isothermal friction coef- 
ficient measured at the same Reynolds number by the following expression 


f'iso 


f 


<1 


\ *^all ] 


(5) 


In figure 5 the nonisothermal friction coefficients fq calculated from 

experimental pressure-drop data are compared with the friction coef- 
ficients fq calculated by means of the above equation. The maximum 

deviation in conparlson is 7 percent, although the deviation for most of 
the data is less than 4 percent. 

Equation (5) is similar to the expressions in the literature which 
apply to water and n-butyl alcohol. The general expression employed in 
the literature for relating fq to fiso 

^iso _ ( ifbu2k\™ 

^q \hrf-all/ 

References 6 and 22 present fluid -friction data on water in which the 
exponent m of the above equation was found to be 0.l4 when applied to 
a range of mean buUt temperatures from 50° to 150° F. Reference 5, which 
is also an investigation with water, indicates that the exponent m is 
0.13 for a range of bulk, tenperatures from 6o° to 125° F. The results 
of investigations with n-butyl alcohol presented in reference 3 show that 
the exponent m is 0.l4 when the ratio of viscosities is less than 3- 
Therefore, the magnitude of the exponent m for WFNA which applies to 
bulk temperatures from 6o° to l45° F is not appreciably different from 
the values reported which apply to water and n-butyl alcohol. The dif- 
ference in the value of m may be attributed to the difference in the 
viscosity properties of the liquids in question. 

The value of the exponent m for WFNA was determined graphically 
from the slope of the logarithmic plot in figure 6. The scatter of the 


^With heat addition. 
^Without heat addition. 




i6 


NACA TN 3181 


data made it necessary to check the maximum deviation of the data from 
a mean value before the slope could be established. The values of the 
isothermal friction coefficients employed in equation (4) were obtained 
experimentally and are presented in figure 7- The curve marked 
"experimental" in figure 7 is a plot of the experimental friction coef'^- 
ficlent of WFNA as a function of the Reynolds number. The other curves 
in figure 7 were obtained from the Moody report (ref. 23). The surface 
finish of the inside of the test section was considered to be similar 
to the finish of commercial drawn tubing. The experimental values of 
the friction coefficient are about 6 percent hi^er than the values 
presented in reference 23 which apply to commercial tubing. 


Correlations of Average Heat-Transfer Coefficients 

Correlation based on average bulk temperature . - The forced- 
convection heat-transfer data for WMA are summarized in figure 8, which 

h 2/3 

presents the Colburn j-factor, j = ^ function of the 

Reynolds number The heat-transfer coefficient h in the afore- 

mentioned relationship for j is the average value for the test section. 
The physical properties of the WFNA used in calculating j and 

were evaluated at the arithmetic average for the bulk temperature. 

Several of the points in figure 8 lie above the curve (straight 
line) recommended Tor correlating the data.- Approximately one-third of 
these test points were obtained when the Inside surface temperature of 
the test section exceeded the saturation temperature of the WFNA so that 
nucleate boiling may have occiirred. Consequently they were not expected 
to correlate with the p-urely convective data represented by the correla- 
tion curve. The points which lie considerably below the correlation 
curve were obtained in run 28 after the acid had been circulated through 
the apparatus for several hours. Figure 9 Resents the deviation of the 
results of runs 27 and 28 from the correlation curve as a function of the 
time that the acid had been circulated through the apparatus. It is seen 
that the heat-transfer coefficient decreased with this time. The reason 
for the significant decrease in the heat-transfer rate is not known 
definitely at this time. Since scale has been reported in boiling heat- 
transfer measurements with WFNA and dissociation reactions at elevated 
temperatures change the properties of the WFNA (ref. 2l), the aforemen- 
tioned deviations from the correlation are attributed to these phenomena. 
The determination of the effect of scale and/or dissociation of the heat- 
transfer rate for WFNA requires further investigation. The data for 
run 28 were extrapolated to zero time and the values corresponding to 
zero time were taken as the j -modulus corresponding to fresh WFNA and a 
clean test section. 
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Figure 10 presents the j -modulus as a function of neglecting 

the runs vhere boiling occurred, and the results of run 28 corrected to 
clean-tube fresh-acid conditions. The data can be correlated by the 
equation 


j = 0.024(%e)-°-^ 


or 


( 6 ) 


( 6 a) 


where all the physical properties of the WFNA are evaluated at its 
average bulk temperature. 

The maximum deviation of the j-iiKDduliis from the correlation curve 
is 32 and -9 percent. The constant in equation (6) can be changed so 
that the maximum deviation would be approximately ±20 percent. Such a 
change, however, places the majority of the points below the correlation 
curve. Equation (6) is, therefore, considered to be the most satis- 
factory correlation equation for the data presented in figure 10. More- 
over, from the point of view of designing a regenerative cooling system 
it gives conservative values. 

The points which lie farthest above the correlation curve in fig- 
ure 10 represent the data taken in nxns 2, 2h, and 29 . Each nn was 
made with a different test section and fresh WFNA and was the first run 
in its respective test section. The test section used for run 24 was 
also used for rims 25 and 26 , but in the latter runs a different charge 
of WFNA was circulated through the apparatus. Since the data for 
runs 25 and 26 correlate satisfactorily and no additional data were 
taken with the test sections used in runs 2 and 29 the discrepancies 
shown by runs 2, 24, and 29 are probably due to differences in WFNA 
composition and not to instrument error. Since WFNA is not a pure 
liquid, but a mixture of nitric acid, nitric oxide, and water, its 
physical properties are functions not only of its temperature but also 
of its conposition. The effect of acid composition on the behavior of 
WFNA should be more thorou^ly and systematically investigated. 

Correlation based on Sleder-Tate equation . - Sieder and Tate (ref. 24) 

proposed that the ratio (m-/ps)^ included in correlation equations 
such as equation (4). This modification is introduced to correct the 
equation for the lack of temperature similarity. Figure 11 presents 
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the par^eter j(Kpe)^*^ as a function of (p/ps) for the data obtained 
in the heat-transfer experiments. Although there is considerable scatter 
in the data, they do indicate that the slope recommended by Sieder and 
Tate, n = 0.l4, is substantially correct. 

Figure 12 presents j/(p/ps)*^‘ as a function of Npg. The best 
stral^t line on the figure for correlating the data appears to be 

= 0.022(NRe)"°'^ (7) 

or 

= 0.022(l^Re)°•®(l^Prf^^(^^^8)°‘^^ (Ta) 


Neglecting the data taken vith dissociated acid or a scaled test 
section (iTin 28), the data correlate within 3^ and -7 percent. Equa- 
tion (7) does not give a better correlation of the experimental results 
than does equation (6). The latter equation has the advantage of 
simplicity. 

Reynolds analogy .- The nonisothermal friction data were correlated 
with the convective heat-transfer data by utilizing the Reynolds analogy. 
En5)loylng the method described in reference 25, which permits applica- 
tion of the Reynolds analogy to liquids, predicted values of the Nusselt 
number were calculated. In figure 13 the experimentally determined 
Nusselt numbers are compared with the predicted values. _The average 
deviation of the calculated values from the experimental values is 
approximately 9 percent. 

In the application of the Reynolds analogy for liquids, one nssump- 
tion was utilized in addition to those employed in reference 25 for 
developing the method. Because of the absence of any data concerning 
the radial temperature distribution of the fluid it was assumed that the 
maximum temperature difference in a radial direction at any station in 
the tube would be the difference between the inside wall temperature and 
the bulk temperature at the station, assuming that the bulk temperature 
varies linearly with the length of the tube. 

Correlation of local values of heat-transfer coefficient .- The 
variation of the heat-transfer coefficient E with length along the 
test section was studied. The 'local values of the heat-transfer coef- 
ficient were measured at seven different locations. At each location 
the physical properties were determined from the local bulk tenperature 
assiunlng a linear variation of the bulk temperature with test -section 
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length. Only those nonboiling trials where the over-all hulk -temperature 
rise for the test section exceeded 15° F were used In order to limit the 
analysis to tests of the higher heat-flux densities. The data from 
rtms 2, 24 , 28, and 29 were omitted because they scattered to such a 
degree that they made no contribution to the analysis. Local values 
of the j -modulus and Resmolds number were calculated for the following 
positions: 


Distance from inlet, 
in. 

Distance/diameter, x/D 

3 

5.6 

6 

11.1 

9 

16. 7 

12 

22.2 

15 

27.8 

18 

33.4 

21 

38.9 


The data were plotted at each location and were correlated satis- 
factorily by the equation 

J = (8) 

where C is a function of the distance from the inlet to the test sec- 
tion and the x/D value. The values of C for the seven locations are 
tabulated below: 


x/D 

C 

5.6 

0.0246 

11.1 

.0241 

16.7 

.0243 

22.2 

.0240 

27.8 

.0234 

33.4 

.0232 

38.9 

.0230 


Figure l 4 presents C as a function of the tube diameter from the 
inlet of the test section. In the same figure are presented results of 
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investigations at Massachusetts Institute of Technology (ref. 6) and 
University of California, Los Angeles, (refs. 7 to 9) on high-pressure 
vater. Since the M.I.T. values of C were hased on physical properties 

determined at the local film temperattire^ tljey are not comparable with 
the values obtained at Purdue University with VTFSA or those obtained at 
U.C.L.A. with water; the latter two investigations based the value of C 
on the local bulk temperature. The trend of the M.I.T. results, however, 
is not influenced by tenperature employed for determining the physical 
property of the water. The U.C.L.A. data indicate that the local heat- 
transfer coefficient decreases rapidly for the first 40 tube diameters 
of the test section. Thereafter it increases gradually. The U.C.L.A. 
investigator attributed the Increase in the heat-transfer coefficient 
after the if-0 diameters to the Increase in liquid tenperature (refs. 7 
to 9). The length of the test section used.ln the WFIfA experiments at 
Purdue University was 44.5 diameters or approximately one-third of the 
x/D ratio of test section used at U.C.L.A. Consequently, the studies 
of the effect of x/d could not be made conpletely conparable with those 
made at U.C.L.A. 

The analysis of the effect of x/D on the heat-transfer coefficients 
indicates that the constant in equation (4) should probably be modified 
to 0.023. Thus 


3 = 0.023(Nf,J 


- 0.2 


(9) 


Equation (9) can then be used to calculate the heat-transfer rate at any 
position, being 'slightly conservative at x/D ratios less than x/D = 50- 


Nucleate Boiling 

Pressure drop with nucleate boiling .- Figure 20 presents the effect 
of heat addition on the pressure drop across the length of the test 
section. In the convective region (without nucleate boiling) increasing 
the rate of heat transmission decreases the pressure drop, but when 
nucleate boiling is encountered the pressure . drop across the test-section 
length is increased. Nucleate boiling appears to begin whenever the 
inside surface temperature of the test section reaches the saturation 
temperature. The region where boiling is assumed to begin is indicated 
in figure l4. No satisfactory analytical method was developed for 
predicting the nonlsothermal friction coefficient in the boiling heat- 
transfer regime. 


^The local film temperature is the arithmetic average of the surface 
and bulk temperature. 
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The pres sure -drop curves for the boiling regime shown in figure l4 
are similar to comparable pressure-drop curves presented in the litera- 
ture (refs. 7 to 10 ). The effect of varying the system pressure was 
not investigated in the boiling regime because the heat -flux density 
required for boiling at pressxires above 100 Ib/sq in. abs could not be 
attained. The data for the curves presented in figure 15 were obtained 
when the system pressure was 6k. J Ib/sq in. abs and the bulk inlet 
temperature was held at appixixlmately 80° F. The upper extent of the 
heat-flux density for each pressure-drop curve was determined by the 
maximum tolerated tube-surface ten^jerature at the downstream end of the 
test section, which was approximately 550° F. 

The rapid rise in the pressure drop across the test section with 
Increasing heat flux after boiling begins is attributed to the occurrence 
of boiling over a greater portion of the inside surface area of the test 
section. The phenomenon of boiling is thought to increase the pressure 
drop because: . 

(a) The nucleate vapor bubbles which form at the inside surface of 
the tube Impose a drag force on the fluid flow 

(b) The formation of vapor at the surface of the tube causes the 
local value of the kinematic viscosity at the wall to increase, signifying 
that the viscous forces of the fluid are becoming more pro min ent at the 
wall 


Nucleate -boiling heat transfer .- Some preliminary investigations 
were made with heat transfer in the nucleate-boiling region. Efforts 
to correlate these heat-transfer results by using methods employed by 
other investigators on different liquids were not successful. 

Several investigators (refs. 5 and J to lO) obtained satisfactory 
correlation of their data in the nucleate-boiling region by plotting 
q/A as a function of the excess temperature Atx = ts - tsat ^ 
constant pressure. To calculate the excess temperature of WFM, the 
vapor-pressure data, obtained at pressures up to 20 Ib/sq in. abs by 
Dr. W. L. Slbbett and others at Purdue University, were extrapolated to 
a pressure of 65 Ib/sq in. abs as shown in figure 16 ; the values of vapor 
pressure based on that extrapolation are open to question. Figure 17 
presents the heat-transfer data where nucleate boiling was indicated; 
q/A is plotted as a function of Atj^. All of the test points on the 
figure except one are taken from run 25. The tests points tabulated in 
figure 17 were not included in the plot because the calculated excess 
temperatures, which were based on the extrapolated values of vapor pres- 
sure, were either zero or negative. It is apparent that no correlation 
of the nucleate -boiling data could be obtained. 
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Another method for presenting nucleate-boiling data is to plot the 
heat-flux density against the temperature difference tg - t-h„Tv for 

constant weight flow rate (refs. 1 and 5 )» it was attempted to correlate 
the nucleate-boiling data for WFNA on that basis. The results are shown 
in figure 18 . The curves dre geometrically similar to those obtained by 
other investigators with fluids other than WFNA. It is seen that for 
each weight flow rate the slope of the Curve increases sharply as the 
surface temperature of the tube reaches the saturation temperature, 
indicating that small Increases in the surface tenqperature result in 
large Increases in heat-transfer rates. The investigations with pure 
liquids show that the data in the boiling region lie on one curve. The 
experiments with WFNA show two curves in the boiling region, the curves 
being approximately 40° F apart. Without more precise data on the vapor- 
pressure-temperature relationship for WFNA, correlation of WFNA data in 
the boiling region cannot be expected. 


Experience with film boiling.- No effort was made to investigate 
the limits of stable heat transfer with nucleate boiling, but film boiling 
was Inadvertently encountered in one trial. With a constant weight flow 
of approximately 2.5 pounds per second, the electrical power supplied to 
the test section was increased to 52 kilowatts which gave an outside 
surface temperatvire of approximately 540° F. The power was then reduced 
to 48 kilowatts which reduced the outside surface temperature to approxi- 
mately 500 ° F. After a very short time the test section failed (bum 
out). The burn-out condition is noted in figure 18 . The tube was burned 


entirely through about 


inches from the outlet end. 

8 


CONCLUSIONS 


The signiflcsint conclusions derived from the study of the heat- 
transfer-and fluid-friction characteristics of white fuming nitric acid 
(WFNA) flowing in an electrically heated horizontal tube are as follows: 

1. Over a range of Reynolds numbers from 55^000 to 220,000 and 
heat-flux densities from 0.13 to 1.4 Btu/(sq in.) (sec) the Fanning fric- 
tion coefficient obtained with heat transfer in the absence of nucleate 
boiling fq can be related to the isothermal Fanning friction coef- 
ficient f^gQ by the expression 


^Iso 



ifbuikr* 

^Vall/ 
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where is the viscosity of the WMA. The average values of heat- 
transfer coefficient for the test section are correlated satisfactorily 
hy the equation: 


J . 0.02l.(Njfe) 


- 0.2 


where j is the Colhurn factor, Nj^ is the Rejmolds ninnher, and the 
physical-property data are evaluated at the buUc temperature of the WFHA. 

2. The local heat -transfer coefficients decreased sli^tly as the 
distance from the inlet to the test section is increased. The decrease 
is attributed to the formation of the thermal boundary layer. For a 
fully developed boundary layer (ratio of distance from inlet to tube 
diameter greater than appMximately 50 ) the following equation gives a 
more satlsfactoiy correlation than does the equation presented under 
conclusion 1 above: 


j = 0.023{Wpe)“°'^ 


where the physical properties of the WFNA are based on the local bulk 
temperature . 

3 . Continued circulation of the same acid in the test apparatus 
reduces the heat-transfer coefficient as much as 25 percent. Whether 
this is due to scale, dissociation, corrosion, or a combination of the 
three is not known. 

k. When the inside surface ten 5 )eratiire of the tube exceeds the 
sat\iration ten^ierature of the WFWA the heat-transfer rate and the pres- 
stire drop are increased over those predicted by the equations presented 
under conclusions 1 and 2 above. No satisfactory correlations of data 
obtained in these preliminary investigations in the nucleate-boiling 
regime were possible; in the case of the heat-transfer data the uncer- 
tainty in the vapor -pres sure data was one of the principal causes contrib- 
uting to the difficulty in correlation. 


Purdue University, 

Lafayette, Ind., February 11, 1953- 
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SPACED ALONE TEST SECTION 
^(2 THERMOCOUPLES PER STATION) 




THERMAL CONDUCTIVITY, BTU/(HR)(FT)(®F) 
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Figure 5 -- Thermal conductivity of Haynes -S tel 1 i 
on data for similar alloys from reference 20 ; 
specific-heat determination . 
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HEAT TRANSFERRED, MEASURED ELECTRICALLY, BTU/SEC 


Figiire L.- Heat balance. 
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EXPERIMENTAL NONISOTHERMAL FRICTION FACTOR 


Figure 5 *- Comparison between calculated and experimental friction factors. 



ISOTHERMAL f 
NON-ISOTHERMALf 


U) 



VISCOSITY at bulk TEMPERATURE 
VISCOSITY at wall TEMPERATURE 


Figure 6.- Relationship betvreen Isothermal f arifi f vlthln 

convectlve-beat-transfer region. Average slope^ 0.12. 
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FRICTION FACTOR 



Figure J.- Isothermal Farming friction factor for a sjzKXjth pipe taken fi*om 

reference 23 . 
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COLBURN j- FACTOR 
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Figure 8.- Correlation of heat -transfer data. Haynes -Stellite Alloy 25 
tube 2k- inches long with 0.559-inch inside diameter; pressure, 64 to 
165 Ib/sq, in. absj t 2 ^, 50° "to 15T® Fj physical properties of acid 
determined at average bulk tengjerature . Line represents equation 

j = 0.024(%e)"°-^- 




DEVIATION FROM CORRELATION j =0.024 (N^g)"®-^ PERCENT 
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Figure 9 »- Effect of scale and/or dissociation on tieat- transfer correlation 

runs 27 and 28. 
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REYNOLDS NUMBER, Nr, 


Figure 10.- Correlation of heat-transfer data. Line represents equation 
j = for conditions, see figure 8. 
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Figure 11.- Influence of viscosity ratio on heat -transfer correlation. 

Slope, O.lli-. 
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REYNOLDS NUMBER 

Figure 12.- Correlation of heat-transfer data using vlscosl-ty-ratio cor- 
rection. Line represents equation - 0'.022 (Hj{q]”®’^; for 

conditions, see figure 8. 
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NUSSELT NUMBER OBTAINED FROM REYNOLDS ANALOGY 
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COEFFICIENT IN COLI3URN EQUATION 



T?*f ~\ li 

J.- W4. ■ 


T /t» 

Xi/U 


xtiuxy oa aea'D-'cransrer correxatiion, JJata 


taken from references 7 to 10 ; coefficient determined "by UCIA for 
KFHA using liquid physical properties determined at hulk temperature 
and hy MIT using film teinperature . 
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Tn^l^'rBTU/(SOIN.)(SEC) 


Figure I 5 .- Effect of heat addition on frictional pressure drop for con- 
vective ftiTd toiling heat transfer, t]^ = 80 ° F; P = 64.7 Ih/sg in. 
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TEMPERATURE DIFFERENCE, 


Figiire I 8 .- Forced convection with and without surface boiling. Haynes - 
Stellite Alloy 25 tube 2k inches long with 0.557-inch inside diameter j 
^bulk^ 91 ° "fco 122° Fj nucleat boiling pressure, 65 Ib/sq in. abs. 
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